Globular cluster luminosity function as distance indicator 
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Globular clusters are among the first objects used to 
establish the distance scale of the Universe. In the 1970- 
ies it has been recognized that the differential magni- 
tude distribution of old globular clusters is very similar 
in different galaxies presenting a peak at My ~ —7.5. 
This peak magnitude of the so-called Globular Clus- 
ter Luminosity Function has been then established as 
a secondary distance indicator. The intrinsic accuracy 
of the method has been estimated to be of the order of 
~ 0.2 mag, competitive with other distance determina- 
tion methods. Lately the study of the Globular Cluster 
Systems has been used more as a tool for galaxy forma- 
tion and evolution, and less so for distance determina- 
tions. Nevertheless, the collection of homogeneous and 
large datasets with the ACS on board HST presented 
new insights on the usefulness of the Globular Cluster 
Luminosity Function as distance indicator. I discuss 
here recent results based on observational and theoreti- 
cal studies, which show that this distance indicator de- 
pends on complex physics of the cluster formation and 
dynamical evolution, and thus can have dependencies 
on Hubble type, environment and dynamical history of 
the host galaxy. While the corrections are often rela- 
tively small, they can amount to important systematic 
differences that make the Globular Cluster Luminosity 
Function a less accurate distance indicator with respect 
to some other standard candles. 
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1 Introduction 

Most astrophysics quantities depend critically on accu- 
rate distance measurements and globular clusters are 
among the first objects used to measure distances. The 
literature on their use as distance indicators is very 
extensive and cannot be reviewed fully here. I opted 
therefore to give a summary on the early history of the 
usage of the globular cluster luminosity function as dis- 
tance indicator (Sec. [2]), before describing the method 
and discussing the way that globular cluster photome- 
try is used to measure distances (Sec. [3]). The main goal 
of this article is to review the recent observational and 
theoretical results and provide an assessment on the 
accuracy and usefulness of the Globular Cluster Lu- 
minosity Function (GCLF) as distance determination 
method. 

Many articles and reviews in the past examined the 
GCLF as distance indicator suggesting the accuracy 
of the method ranging between 0.1-0.2 mag, at the 
level comparable to, or sometimes even bette r than 
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dicator bv lRichtlerl |2003l) examines in particular those 
galaxies that yield discrepant distance measurements 
when the GCLF distances are compared to those of 
surface brightness fluctuation method, concluding that 
discrepant cases are most probably caused by peculiar 
globular cluster populations that are not uniformly old, 
and therefore inapplicable for the method. The most 
recent critical comparison of the relative Virgo vs. For- 
nax cluster distance as measured using GCLF and other 
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Fig. 1 The Milky Way (top) and M31 (bottom) globu- 
lar cluster luminosity fu nctions based on data published by 
iKron and Mavalll l|l96Cl ) in Tables XV and XX, respectively. 
In the sample of the Milky Way globular clusters, the abso- 
lute magnitudes for th e two brightest clust e rs (47 Tuc and 
ui Cen) are taken from lGascoigne and Burr] (Il956l), as men- 
tioned in a note added in the proofs by Kron and Mavall 
l|l96Cl ) 



d2010l) 



light, which they interpreted as the probable range of 
distance moduli for M31. More importantly, the au- 
thors concluded that "the present data do not indi- 
cate any substantial difference ... between M31 and 
the Galaxy." The luminosity functions of the MW and 
M31 clusters based on data published in that work are 
shown in Figureffl The MW GCL F is drawn from Table 
XV from lKron and Ma vall (I960 ), augm e nted w ith the 
measurements from Gascoigne and Burr ( 1956| ) for the 
two most luminous clusters, 47 Tuc and ui Cen (the un- 
intentional omission of these two clu sters is mentioned 
in a note bv lKron and Mavalll (|l96fj ). The M31 GCLF 



standard candles is presented by IVillegas et al 
and is discussed briefly here. 



2 A historical note 

Old globular clusters are ubiquitous around all large 
galaxies. They are compact stellar systems containing 
10 3 to 10 6 M Q in stars, and spanning a range of mag- 
nitudes from My ~ 2 to as bright as My = — 11 mag. 
Since they are among the highest luminosity objects in 
galaxies, they have been recognized in the early 1960-ies 
as a potentially important distance indicator, and even 
before that in the second decade of the 20th century 
Shapley used them to investigate the size of t he Milky 
Way gal axy (MW) and our position within it (|Shaplevl 
1918alibl ). 

In their 1960 paper on "Photoelectri c photometry of 
galact ic and extragalactic star clusters" Kron and Mav allj 
(1960) derived the photographic and visual magnitude 
distributions of globular clusters in the Galaxy and 



M31, adopting the phot ometric magnitude M pn = 0.0 
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for RR Lyrae variables ( Sawyer Hogglll959h . The two 
distributions were found to have a constant difference 
of 23.6 mag in photographic light and 24.1 in visual 



is based on their Table XX. The absolute magnitudes 
are d erived assumin g a distance modulus of 24.6 for 
M31 (ISandagelll958l) . The hatched histograms show the 
distribution for outlying clusters, which are presumably 
less affected by rather uncert ain extinction corrections . 
Based on these observations Kron and Mavall (1960) 
derived the mean magnitude of the globular clusters in 
the Galaxy of (M) Pt = -7.7 and (M) vt = -8.2 mag. 
Applying these values to the observed mean magnitude 
for 11 clusters in the Large Magellanic Cloud their esti- 
mated distance modulus for the LMC was 19±0.5 mag. 

In an indep endent study of M31 globular clusters 
Vetesnik ( 19621 ) noticed the similarity of M31 and MW 



globular cluster luminosity distributions with a peak at 
My ~ —8 mag. It should be noticed though that these 
early studies suffer from significant biases due to incom- 
pleteness and primary distance calibration uncertain- 
ties. In both galaxies the observed clusters are almost 
all brighter than My ~ —7. In addition, the corrections 
for often unknown extinction are not always applied. 

Beyond the Local Group one of the first galaxies to 
have its distance determined using the observations of 
its globular clusters is M87. Two studies published in 
1968 derived values for the Hubble constant based on 
the derived distance to M87, in remarkable agreement 
with mo dern measurements ma de within t he HS T Key 
Project ( Freedman et al. 2001 ). Racine ( 19681 ) mea- 



sured magnitudes of more than 1000 globular cluster 
candidates around M87 and derived a distance modu- 
lus of (m - M) B = 30.7 ± 0.2 to M87 based on the 
comparison of the brightest cluster magnitudes with 
the brightest clusters in the MW and M31, as well as 
based on the peak magnitude of the cluster distribu- 
tion. His value of the Hubble constant (he actually calls 
it "a velocity /distance ratio") is 79 ± 12 km/sec/Mpc. 
In a similar study entitled "A New Determination of 
the Hubble Constant fro m Globular Clu sters in M87" , 
published the same year, Sandagel ( 1968 ) derived a dis- 
tance modulus to M87 of (m — M)ab — 31.1, based on 
brightest cluster magnitude, and the Hubble constant 
of 75.3±^ km/sec/Mpc. 



3 



While these early studies suffered from limitations, 
the prominent one being the inconsistency whether to 
use the brightest cluster magnitude, the average mag- 
nitude of a sample of globular clusters, or to derive the 
peak magnitude in the distribution, they have already 
in 1960-ies demonstrated that globular clusters may be 
useful distance indicators. Some more modern deter- 
minations of the Hubble constant using the GCLF are 
discussed in Section [6] 

In 1970 Racine remarked that the brightest clusters 
in M87 and in the Loc al Group gal axies, the MW and 
M31, may be different (|RacineHl97olh and shortly there- 
after the brightest clusters were not used any more as 
standard candles. Instead the GCLF maximum value 
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was established as the standard candle (|Haneslll977aj) 
and this method is described in the next section. 

Given that the brightest clusters are not a standard 
candle, I do not discuss them further here. It is worth 
mentioning, however, that they are well worth a study 
on their own, because they contain important clues for 
the GC LF formation and overall shape ( see th e discus- 
sions bv lMcLaughlinlll994tlRTchtlerl2003l[2006l) . and be- 



cause they represent a link between the star clusters and 



Transition Obiects (IHilker et al. 1999; Haseean et al. 


2005; 


Reikuba et al. 12007; 


Mieske et al. 12002. 20081). 



3 The GCLF distance determination method 

The main premise of the GCLF method is that the 
old globular clusters magnitude distribution exhibits a 
universal shape that is characterised by some number 
that can be used as the standard candle to measure 
the distances. The distance determination method is 
then conceptually very simple and observationally rel- 
atively inexpensive, as it does not require time-series 
imaging: it is only necessary to measure luminosities 
of all globular clusters in a given target galaxy reach- 
ing a sufficiently deep limiting magnitude to recover the 
characteristic magnitude, build its magnitude distribu- 
tion, and compare it with the known magnitude dis- 
tribution of globular clusters in the MW for which the 
GCLF standard candle is calibrated based on accurate 
individual cluster distances from well understood and 
well calibrated primary distance indicators. Of course, 
the real measurement has additional complications, and 
several corrections and sources of error are discussed 
below and in Section [7] 

Figure [5] shows the modern version of the GCLF for 
the MW globular clusters based on the 2010 compila- 
tion of the Ha rris McMast er Milky Way Globular Clus- 
ter catalogue (|Harrislll996h . The black histogram shows 
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Fig. 2 T he Milky Way GCLF based on the 2010 ver- 
sion of the iHarrisl l| 19961 ) catalogue for 129 clusters with 
E(B — V) < 1.0 is shown with black histogram. The subset 
of these clusters that have [Fe/H] < — 1 dex is plotted as 
hatched blue histogram 



the distribution of 129 clusters with E(B - V) < 1.0, 
while a subset of these that are more metal-poor than 
[Fe/H] < — 1 dex is shown as hatched histogram. Over- 
plotted to both is a Gaussian fit for the distrib utions. 
As was already men tioned by other authors (e.g. lHarris 
200ll; iRichtlerlliooj ) the name Globular Cluster Lumi- 
nosity Function is a misnomer due to the fact that it 
refers to the magnitude and not luminosity distribution 
of globular clusters in a galaxy. 

The most commonly adopted way to derive the char- 
acteristic magnitude is to fi t the GCLF wi th a Gaussian 
function, as introduced by lHanesI (jl977bh : 



dN 



1 



dM v aVZn 



exp 



(1) 



The distribution is then characterized by only two pa- 
rameters: the dispersion (sigma) and the peak (or 
turnover, TO) magnitude. The latter is the standard 
candle and the distance measurement is done relative 
to the MW or to the M31 GCLF TO value. The GCLF 
is a secondary distance indicator, since the absolute dis- 
tances to the MW and M31 globular clusters must be 
known. These are based primarily on the RR Lyrae 
luminosity scale for the MW globular clusters, or alter- 
natively on the Cepheid distance to M31 (Section [5]). 

It has to be emphasised that the fitting of the glob- 
ular cluster magnitude distribution with a Gaussian 
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Fig. 3 The Milky Way GCLF in magnitude units ( left) and in luminos ity units (right) fitted with the Gaussian (solid 
line) and t$ function (dashed line). This figure is from lLarsen et al l l|200ll) Figure 1, reproduced by permission of the A AS 



provides just a useful parameterisation with a n ana- 
lytic function, but it is not physically motivated ([Harris! 
2001). This shape of the GCLF is a consequence of the 
logarithmic magnitude scale combined with a change 
in power law slope of the underlying distribution of lu- 
minosities of globular clusters, which, combined with 
the approximately constant mass-to-light ratio for old 
globular clusters, leads to a power-law mass spectrum 
of the globular clusters in the form N(M) oc M~ a 
( Harris and Pudritj 1994 ). At the magnitude that cor- 
responds to the peak of the Gaussian the linear lumi- 
nosity function (dN/dL oc L~ a ) changes the power-law 
slope from a < 1 to a > 1, because then the magnitude 
distribution exhibit s t he peak dN/ d rriy oc io 4 ( 1 - Q ) m ^ 



(jMcLaughlinl 11994) . llVlcLaughlinl (|1994l) showed that 
actually the GCLF is non-Gaussian, but as it is ob- 
served to be peaked, unimodal and nearly symmetric, 
its approximation with a Gaussian provides relatively 
unbiased distance estimate. 

Alternative GCLF fitting functional forms for the 
MW and M31 cluster distributions include for exam- 



and a piecewise ex ponent i al fit ([McLaughlin 



Baum et al.l 119951 ) . ISeckerl (|1992l ) used a maximum 
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likelihood estimate to demonstrate that the best fit to 
the MW and M31 GCLFs is obtained with the t 5 dis- 
tribution in the following form: 



dN 



dm SV^crt 



1 



(to — too)'' 



(2) 



which presents a longer tail towards the faint end. The 
more extended faint end of the GCLF is visible in the 
MW GCLF (Fig. [5]) and is even be tter appreciat e d from 
Figure[3l which is reproduced from lLarsen et al.l (|2001 ). 
The deviation from the symmetry at the faint end of the 
GCLF is a consequence of the dynamical evolution of 
the initial cluster luminosity function which is expected 
to be better represente d by a power -law than a log- 
normal distribution (e.g. lLarsenll2002l) . In spite of this 
difference at the faint end, the fit of the GCLF TO 
yields very similar value irrespective of w hether the fit 
is don e using a Gaussian or a t$ function (jLarsen et al 
200ll) . 



pie Cauchy and t distribution ( Seeker 1992) , Gauss- 
Hermite expansions (|Abraham and van den Berghll9 95)| 



Figure 0] (adopted from iJordan et al.l (|2007l) ) shows 
the physically motivated functional form for the evolved 
Schechter mass functions dN / dM and the correspond - 
ing GCLFs dN/dm as proposed bv ljordan et al. ( 2007 ). 
The mass distribution of young star clusters in a galaxy 
is well described by a power law function dN /dM oc 
M~P with a characteristic slope of (3 ~ —2. At the high 
mass end a truncation of the star cluster initial mass 



function is observed (e.g. M 51: Gieles etlrt"1l2006al) . (the 



Antennae: I Whitmore et al.lll999[ h" This can be well de 



Burkert and Smith] (j2000h . In a functional 



scribed with an exponential cut-off similar to what was 
proposed for the initial lumino s ity fun ction of globular 
clusters by 
form of the 

index of ~ 2 and an exponentia l cutoff above some large 
mass scale (jGieles et al.ll2006bl ) the initial cluster mass 
function is: 



Schechterl (|19761 ) function with a power law 



dN 



Kissler et al.l (|1994) also prefered t he function a l form — — oc M Q exp 



m 



for the GCLF of NGC 4636, and lLarsen et al.1 (|2001 ) 



dM 



(3) 



used it to fit GCLFs of 17 nearby galaxies. 

The superior fit by the t 5 distribution with respect 
to the Gaussian is due to the asymmetry of the GCLF 



Due to dynamical erosion of preferentially low mass 
clusters the original rising low mass end of the young 
clusters distribution turns down below the power-law 
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Fig. 4 The evolved Schechter mass functions dN/dM 
are shown in the top panel and the corresponding GCLFs 
dN/dm are in the bottom panel. This is Figure 1 from 
I Jordan et al.l |2007t ). reproduced by permission of the A AS 



behaviour for clusters with masses smaller than about 
2 x 10 5 M , which is the mass of a typical globular clus- 
ter at the TO magnitude. Such an evolved Schechter 
mass function has the form: 



dN 
dM ^ 



A) 



"exp 



(4) 



where the cumulative mass loss rate is A = fi ev t such 
that each cluster has a time dependent mass M(t) = 
Mq — 5. Relating this function to the GCLF, by assum- 
ing a typical mass-to-light ratio for globular clusters, 
the model GCLF has the following form: 



dN 



-j^Q— 0.4(m— m c ) 



dm (io-o.4(m-m c ) + io-°- 4 («-"*<0) 2 



exp 



(-10" -' 



where the magnitude m = C — 2.5 log m, 5 = C 
2.5 log A, and M c = C - 2.5 log M c . 



Harris et al.l f|2009h note that there is no distinguish- 
able difference between Gaussian and Schechter-likc 
function fits for the measurement of the GCLF TO 
magnitude. 



4 Theoretical models 

Before discussing the calibration of the GCLF TO mag- 
nitude, I provide a very brief overview of the disruption 
processes that erode the initial globular cluster mass 
function (GCMF) to the present day GCMF that is re- 
lated to the observed GCLF. 

The dynamical processes that destroy globular clus- 
ters in a fixed galactic potential are (i) mass loss due 
to stellar evolution (this process does not change the 
shape of the mass function), (ii) dynamical friction 
(Tdf oc M _1 ), (iii) tidal shocks heating by passages 
through bulge or disk (r s h oc phPcr where ph is the 
mean density within the half mass radius Rh), and (iv) 
evaporation due to 2-body relaxation (r ev cx M/^fph). 
In particular the last two are most important for the 
transformation of a power la w or Schechter-like initial 
GCMF described by Eq. GH (|Harris and Pudrit j Il994 : 



Gieles et a l. 2006a) 



The p io neerin g studies of Fall and Zhana (I200lh 



Vesperinil (|200l[ ): iMcLaughlin and Falll (|2008[ ) suc- 
cessfully reproduced the shape of the GCMF by 
cluster evaporation. However, the fundamental as- 
sumption in those studies is the mass-independent 
mass-loss-rate, which has been shown to be incor- 



rect ( Baumgardt and Makino 2003; Lamers et al. 


2005; 


Kruiissen and Porteeies Zwart 


2009 


). In particular, 
) demonstrate that 


Kruiissen and Porteeies Zwart 


(2009 



(5) 



due to preferential depletion of low mass stars from 
globular clusters, the M/L ratio of the clusters changes. 
They model the influence of the luminosity dependency 
of M/L on the globular cluster mass function derived 
from the obse rved GCLF, providing an improvement to 
the model of Fall and ZhansJ (| 2 lh , for the first time 
fitting the complete GCLF (and the derived GCMF) 
with a model (Figure [5]). 

An alternative view of the evolution of the glob- 
ular cluster mass function and evidences for an ap- 
proximately Gaussian initial mass function of old glob- 
ular c lusters is presented by Parmentier and Gilmorel 
( 20071 ). who stress the importance of the gas removal 
process and the protoglobular cloud mass scale. 

The next step is the full blown modeling of the 
star cluster populations formation and evolution within 
galaxy simulations, that have changing potential due 
to their merging his tory. The first results from 



Kruiissen et al.l (|201ll ) provide a better understanding 



of the cluster evolution and their observed properties. 
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and evolution histories, and as a possible observational 
evi dence of the hierarchica l assembly of lar ge galax- 



log(L v /L ) 



Fig. 5 The comparison of the observed (filled) 
and modeled (dotted) GCLF. Th is is Figure 3 from 
iKruiissen and Portegies Zwarj ||2009h . reproduced by per- 
mission of the AAS 



5 The calibration 

Irrespective whether the GCLF TO is fitted with a 
Gaussian, t$ , or an evolved Schechter function, the fun- 
damental question that needs to be answered before it 
can be used as a standard candle is whether the GCLF 
has the same TO peak value in different galaxies. This 
has been addressed in numerous papers since the early 
days of the application of the method. Here only few 
selected recent papers are mentioned in which the TO 
value was examined and calibrated. 

In their 1992 critical review of select ed techniques 



for ex tragalactic distance measurements iJacobv et al 



ies jAshman and ZepJll993 ICote et al.lll998h . There 



fore from the early 90-ies many studies of the globular 
cluster systems in external galaxies were concentrat- 
ing on deriving the color, metallicity, kinematics, spe- 
cific frequency, sizes and age properties for the two 
population s, rather tha n employing them as distance 
indicators ( Harris! 2001 ). The presence of two globu- 
lar cluster populations may have brought into question 
the possibility to obtain a unimodal universal globu- 
lar cluster luminosity function, and therefore relatively 
less effort was put in accurate and precise calibrations 
of the GCLF as dist ance indicator. 

Whitmore] (1997) reviewed globular clusters as dis- 



(|1992| ) concluded that the GCLF method appears to be 
more accurate than expected, but possibly has a scale 
error in the sense of overestimating distances by about 
13%. That paper emphasised the need to collect high 
quality observational data for a list of selected large 
disk and elliptical galaxies close enough to clearly mea- 
sure their GCLF peak magnitudes in order to provide 
a robust calibration of the standard candle, as well as 
the need for a theory of cluster formation and dynam- 
ical evolution that can predict the full distribution of 
globular cluster masses and luminosities. 

Ro ughly at the same time of the review of Jacobv et al.^ 
(1992), in the beginning of the 90-ies, a remarkable new 
discovery in the field of extragalactic globular cluster 
systems studies was the bimodal co lor distribution of 
globular clusters in giant galaxies (IZepf and As hma: 



tance indicators concluding that the intrinsic disper- 
sion in Mv.to is ~ 0.12 mag for bright ellipticals 
and providing the mean value of the Gaussian dis- 
tribution (M v>TO ) = -7.21 ± 0.26 with the width 
a = 1.35 ± 0.05 mag based on literature compilation 
of GCLF measurements. 

Anot her extensive discussion of the GCLF calibra- 
tion by Harris! ( 2001 ). again based on literature se- 
lection of reliable GCLF measurements, presented ev- 
idence for a small difference between the average TO 
magnitude for giant ellipticals and disk galaxies. The 
sample of 16 giant elliptical galaxies has (Mv,to) = 
—7.33 ± 0.04 and mean a of the Gaussian fit of (a) = 
1.4±0.05, while in 14 disk galaxies the mean TO param- 
eters are (My,xo) = -7.46 ±0.08 and (a) = 1.2 ±0.05. 
The accuracy of this calibration, however, is limited by 
the accuracy of the distances to SO and elliptical galax- 
ies in Virgo, Fornax, and Coma clusters, Leo group and 
other nearby groups, as well as the accuracy of the pri- 
mary distance indicators RR Lyrae and HB luminosity 
for the MW, and Cepheids and possibly the tip of the 
red giant branch distan ces to M31. 



1993t I Ashman and" ZepJ 19981: iBrodie and Straderll200 
Kissler-Patia |2000|) . The presence of the two episodes 
of vigorous star formation evidenced by the bimodal 
globular cluster population was quickly recognized as 
the valuable tool to investigate giant galaxy formation 



In the same review lHarrisI (|200l[ ) examines in de- 
tail the MW GCLF, and provides, after comparing the 
MW and M31 GCLFs the TO magnitude of My = 
—7.68 ± 0.14 for the MW projected halo sample that is 
constructed to best represent the globular clusters sam- 
ple as it would be observed in other galaxies. The TO 
magnitude for the M31 GCLF halo sample reported in 
the same review is M v = —7. 80 ±0.12. While the val- 
ues for these two closest galaxies are consistent within 
their respective errors, they appear to be significantly 
brighter than the average for 14 disk galaxies. Further- 
more, the fact that the single review, which admittedly 
was not focused on the GCLF as the distance indicator 
primarily, though it did discuss it at a certain length, 
provided several different values for the TO magnitude 
of the GCLF, later led to different authors adopting 
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Fig. 6 Comparison of the cumulative luminosity function 
for NGC 5128 (thick red line) with that of the Milky Way 
(thin green line) and the M31 (dotted blue line) for V (left ) 
and U band (right). This is Figure 9 from IReikubal |200ll ). 
reproduced with permission A&A ©ESO 



different calibrations in their works. This and other in- 
dependent calibrations resulted sometimes in different 
distance measurements based on the same observational 
datasets. 

In addition to the calibration of the TO the ques- 
tion of the universal shape of the whole GCLF was 
addressed in a number of papers. Most of them com- 
pared in some detail the MW and M31 GCLFs and 
this exercise was repeated as new, d eeper data be- 



came available (typically d ata fo r M31; van den Berghl 



19851: iRacine and Harris] Il992t 



Barmbv et al.1 2001; H uxor et al 



Ashman et al.l 



199 



5; 



20111) . While the re- 



sult depends on the quality and depth of the photome- 
try, and the sample selection, the Kolmogorov-Smirnov 
tests show that the M31 GCLF is slightly brighter and 
narrower than that of the MW, but the differences are 
typically not statistically significant, and can be ex- 
plained by a different metallicity of the two popula- 
tions. 

However, the most common application of the GCLF 
method is to early type galaxies. Elliptical and SO 
galaxies are particularly suitable targets because of the 
large numbers of old globular clusters in their halos. 
Additionally, in spiral galaxies contamination from the 
young star cluster populations, as well as crowding and 
extinction correction problems, are much more severe. 
Hence for a long time the GCLF distance determina- 
tion method relied on a leap of faith that the luminos- 
ity functions of old globular clusters did not depend 
strongly on the Hubble type. 



( ReikubaluOOll see also Fig. [6]). In spite of the signif- 
icant efforts invested, the very large extension on the 
sky combined with the need to spectroscopically con- 
firm the cluster candidates due to significant contami- 
nation from the MW foreground stars in typical ground 
based searches, NGC 5128 still does not have a deep 
and complete enough GCLF suitable for an accurate 
distance determination. Apart from using the globular 
clusters as galaxy formation and evolution probe, the 
GCLF in this galaxy may be of particular interest as 
it could in principle provide an important independent 
cross-check of many different distance dctcrminators. 
NGC 5128 has dista nces determined fr om Cepheids 



( Ferrareseet aJ J 120071). tip of the RG B ([Harris et al 



1999t iReikuballiool lHarris et al1l2010l). p l anetary neb 
ulae luminosity function (|Hui et ah 19931 lHarris et al 



I2010I) . Mira period-luminosity relation (IReikubal 1200 



and s urface brightness fluctuation method (jTonrv et al 
200 ll ). In addition novae have been observed in this 



galaxy (jCiardullo et al.lll990l) and it has even hosted a 
quite typical type la SN (jCristiani et al . 1992). Extend- 
ing the baseline of targets beyond the LMC and M31, in 
which cross-checks of different distance determination 
m ethods can be direct l y com pared, is important. 



Di Criscienzo et al.l ([20061 ) used several recent cali 



brations of the absolute magnitude of RR Lyrae stars 
(My (RR)) to calibrate the TO V-band magnitude of 
the MW GCLF. The derived peak magnitude depends 
strongly on the selection criteria for the MW glob- 
ular cluster sample, as well as, to a lesser extent, 
on the adopted My (RR)vs.[Fe/H] relation. Draw- 
ing the full cluster sample from the H arris McMas - 
ter Milky Way glob u lar clu ster catalogue (|Harrislll996l ) 
Di C riscienzo et al. d2006l) derive the peak magnitude 



My (TO) = -7.40 ± 0.09 mag. This va lue is identica l 
to the one often quoted as coming from Harris! ( 2001 ) . 
It should be noticed though that the lHarris! (|200ll ) re- 
view shows the Gaussian curve with M v = —7.40 and 
a = 1.15 in Figure 1.33, prior to the detailed dis- 
cussion of the GCLF as distance indicator. In the 
later sections this review derives My = —7.68 ± 0.14 
for the MW projected halo sample that is constructed 
to best represent the globular clusters sample as it 
would be observed in other galaxies. This is very 
similar to My (TO) = - 7.66 ± 0.11 mag derived by 
Di Criscienzo et al. ( 20061 ) for a sub-sample of MW 



Th e nearest early type galaxy is NGC 5128 ([Harris et 
20ld) . and its GCLF has been sampled to 2.5 mag 
fainter than the TO (but only in two fields covering a 
small fraction of the total spatial extent of the galaxy) , 
allowing to demonstrate that the difference between 
the GCLF of NGC 5128 and the MW is not larger than 
the difference between the M 31 and the MW GCLFs 



L£ 

al.( 



lobular clusters with E(B - V) < 1.0 and 2 < R GC < 
5 kpc. The difference between the TO magnitudes of 



the two samples is due to differences in their average 
metallicity (the latter being more metal poor), as well 
as to environmental effects, since the two cluster groups 
have a different dynamical history (having been drawn 
from different parts of the Galaxy). 
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Fig. 7 The shift in the peak of the GCLF as a function of 
metallicity for B,V,R,I and J filters. This figure is based 
on sim ulation results reported in Table 3 of I Ashman et al] 



The shift in the peak of the GCLF as a function 
of metallicity for the B, V, R, J, and J filters has been 
investigated bv lAshman et al.l (|l995l) based on simula- 
tions that included the cluster masses drawn from a par- 
ent population similar to that of the MW, and a mean 
metallicity ranging from —0.2 dex to —1.6 dex. The 
GCLF TO in the V-band is brighter for more metal- 
poor globular cluster systems, while in near-IR filters 
the more metal-poor systems have fainter GCLF TOs. 
The smallest difference is reported for the /-band where 
a 1.4 dex shift in [Fe/H] results in only 0.17 mag differ- 
ence in the peak ma g nitude . Figure [7] based on Table 3 



from Ashman et al.l ( 19951 ) summarizes the results of 



these simulations. It shows the GCLF TO magnitude 
shift as a function of metallicity in agreement with the 
observations. Applying the selectio n to the metal-poor 



and m etal- rich sample in the MW iDi Criscienzo et al 



(2006) find the metal-poor globular clusters sample 

([Fe/H] < -1.0, ([Fe/H]) 1.6 dex) systematically 

brighter by V = 0.36 mag with respect to the metal-rich 
sample ([Fe/H] > -1.0, ([Fc/H]) 0.6 dex), and the 



same difference is found for M31 clusters ([Barmbv et al 
200lt 



Di Criscienzo et al. ( 20061 ) adopted the same dis- 



tance modulus to the LMC of 18.50 mag for both the 
RR Lyrae and the Cepheids distance scales, and used 
the latter to set the absolute magnitude of the surface 
brightness fluctuation distance scale. Hence, based on 



the same LMC fiducial distance scale they compare the 
TO magnitudes of metal-poor cluster luminosity func- 
tions in the MW, M31 and the more distant galaxy sam- 
ple of lLarsen et al.1 (|200lh . They find them to be in per- 
fect mutual agreement with My (TO) = -7.66 ± 0.11, 
—7.65 ± 0.19 and —7.67 ± 0.23 mag, respectively. 



6 Potential for measurement of the Hubble 
Constant 

H was derived using the brightest globular clusters as 
soon the individual cluster mag nitudes were measured 
aroun d M87 in late the 60-ie s dRacine! 1 19681 : ISandagel 
19681) . Ivan den Bergh et al.l jl985j) used the GCLF 



TO magnitude of M87 to obtain the Hubble param- 
eter value of 68 ± 10 km s _1 Mpc -1 , although in that 
work the authors emphasise that the assumption of the 
invariant GCLF is "as yet untested" . 

Often the main goal of the globular cluster system 
studies was the investigation of the parent galaxy for- 
mation and evolution history, and GCLF distances were 
done as a secondary use of the data. The uncertainty 
sources are discussed i n the next section, bu t just to an- 
ticipate as stressed by iKissler-Patisd (|2000h : in partic- 
ular background/foreground contamination, aperture 
corrections, incompleteness corrections, GCLF fitting 
errors, possible contamination from a second younger 
cluster population, and uncertainty in the absolute cal- 
ibration of the GCLF TO, were not always fully taken 
into account, re sulting in large scatter and inhomoge- 
neous distances (jFerrarese et al.ll2000al ) . Therefore, the 
compilation of all Hubble constant measurements using 
the GCLF method is beyond the scope of this review. 
The GCLF Hubble diagrams based on selected sampl e 
of galaxies a re sho wn in reviews b y Whitmorel (1997), 
Harris! (|200ll) . and iRichtlerl j2003l) . 



Until today, the universality of the GCLF is de- 
bated and de b atable (see in particular the work of 
IVillegas et al. ( 2010l ) discussed in the next section). 
With deep HST imaging, reaching V ~ 28 or fainter, 
the GCLF method, due to the brightness of its 
TO magnitude, has the potential to step directly 
into the Hubble flow. Hence, it can provide mea- 
surements of the Hubble constant, without passing 
through the next distance ladder rung. It can also 
be used t o verify SN la distances independent ly from 
Cepheids (jdella Valle et al.lll998l : lRichtlerll2003l ). While 
most Hubble constant derivations from the GCLF 
method result in the range of ~ 62 — 83 km/s /Mpc 



(Whitmore et al. 1995; Sandage and 



Whitmorel li997UHarrisll200ll; IRichtlerl 120031) . thev 



ammann 1996 
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Fig. 8 2-band GCLF turnover magnitude vs. z-band 
galaxy absolute magnitude for all galaxies in the ACS Virgo 
(black symbols) and ACS Fornax (red symbols) cluster sur- 
veys. The lines correspond to simultaneous error-weighted 
linear fits performed over both samples including an offset 
of A(m — M) = 0.20 ± 0.04 mag, which is the relative dis- 
tanc e modulus between t he two clusters. This is Figure 6 
from lVillegas et all (|2010T l. reproduced by permission of the 
AAS 

primarily based on Virgo cluster (and M87) GCLF dis- 
tances that st ill require significant co rrect ion for reces- 
sion v elocity. iKavelaars et al. (120001) and lHarris et al 
(2009) measured GCLF distances to Coma giant el- 
lipticals obtaining Hq of 69 ± 9 and 73 km/s/Mpc 
respectively. 



7 Uncertainties 

Without a discussion of uncertainty sources, this review 
would be incomplete, though I point out t hat this topic 
has be e n covered very w el l by rev i ews o f iJacobv et al. 
( 19921 ); IWhitmorel dl997h ; Irlarrisl f|200lh and most re- 
cently bv lRichtlerl (|2003l ). 

Based on the most recent, thorough and homoge- 
neous study, that used highest quality deep HST ob- 
servations of 89 early-type galaxies in th e Virgo cluster 
and 4 3 galaxies in the Fornax cluster, IVillegas et al. 
( 20ld ) derive the relative distance modulus between 
these two galaxy clusters A(m — M) = 0.20 ± 0.04 mag 
(see also Figure [5] adopted from that paper). This value 
is systematically 0.22 mag lower with respect to rel- 
ative distance modulus obtained from surface bright- 
ness fluctuation measurement on the same dataset 



(jBlakeslee et al.ll2009[) . Comparison with the relative 
distance modulus obtained with the planetary nebu- 
lae luminosity function (PNLF) method shows also a 
systematic offset of 0.1-0.15 mag, and the offset with 
respect to Cepheids is 0.27 mag, again for both in the 
sense that the relative GCLF distance Virgo- Fornax is 



small er than for PNLF and Cepheids (|Villegas et al 
2010). This points to a systematic error inherent to 



the GCLF method that makes this distance indicator 
less accurate with respect to the SBF, PNLF or Cepheid 
distance indicators. 

All possible sources of error are considered here, di- 
vided into two groups. 

7.1 General uncertainty sources 

With high precision, deep CCD photometry, the er- 
rors inherent to any photometric method amount to 
~ 0.05 mag and include the following: 

1. Photometric zero point calibration and photometric 
magnitude uncertainties for the observed clusters. 

This error depends on the telescope, instrument, 
and the filter system used. With well understood 
large telescopes and in particular the stability of 
the HST, accuracy at the percent level is achievable. 
Here one should include the aperture correction un- 
certainty for slightly resolved clusters. 

2. Uncertainty associated with extinction correction. 

The GCLF method is principally applied to giant 
elliptical galaxies which have little or no internal ex- 
tinction, and thus the only important contribution 
is from the foreground extinction correction error. 

3. Systematic errors in the fitting procedure and deter- 
mination of the GCLF TO magnitude (mo). 

While the different forms of the fitting functions 
were discussed above, the fitting method, e.g. direct 
Y 2 fitting or maximum li kelihood, can bias the result 



(|Secker and Harris! Il993h . Other issues to consider 
here are the photometric limit and completeness cor- 
rections, as well as corrections for the foreground 
contamination by stars and background contamina- 
tion from compact galaxies masquerading as glob- 
ular clusters. In galaxies within ~ 5 Mpc, under 
best sub-arcsec seeing conditions, it is possible to 
distinguish globular clusters from f oreground star s 
due to their slightly resolved profiles ( ReikubalboOl ) , 
while HST can resolve cluster s up to about 40 Mpc 
( Kundu and Whitmorell2001allbl ; Ijordan et all 12009 : 
Harris) 2009 ). More distant galaxies need statisti- 
cal subtraction of the contaminants. The size and 
quality of this correction dep ends on the filter sys- 
tem used (jDirsch et al.ll2003l ). on the (Galactic) co- 
ordinates of the target, and the availability of the 
suitable, equally deep, comparison background field. 
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7.2 GCLF specific uncertainty sources 

The errors particular to the GCLF method amount to 
~ 0.3 mag and include the following: 

1. Uncertainty in the primary calibrator luminosity - 
distance scale. 

The individual globular clusters in the MW, 
the primary GCLF calibrator, are based on pri- 
mary distance indicators such as RR Lyrae, hori- 
zontal branch, main sequence, sub-dwarf, or white 



(e.g. 


Roes 


19961 Gratton et alJ 


1999: Zoccali et al. 


2001: Richer et all 12004 Recio-Blanco et al. 2005; 


Coppola et al. 


l201ll). For M31, an independent 



GCLF calibration is possible using Cepheids and 
the tip of the red giant branch. Their accuracy is 
thought to be at the level of ~ 10% (for details I 
refer to other reviews in this volume). 
Intrinsic dispersion of mo and the dependence on the 
gl obular cluster sample . 

Villegas et al. ( 2010l ) demonstrate that the intrin- 



sic dispersion of mo is 0.21 mag for the Virgo cluster 
and 0.15 mag for Fornax. 
Environment and dynamical evolution. 

This is the domin ant uncertainty of the method. 
Villegas et all (|201dh find a systematically lower rel- 



ative distance between Fornax and Virgo clusters 
using GCLF when compared to SBF, PNLF, and 
Cepheid distances. This may imply that galaxy clus- 
ters with higher velocity dispersion (higher masses) 
host galaxies with fainter TOs i n their globular clus- 
ter sy stems as suggested by Blakeslee and Toiirvl 



(199 



mo dependence on the Hubble type or galaxy lumi- 
no sity. 

Harris! ( 2001 ') finds a dependence on the Hubble 



type a t the level of ~ 0.13 mag, while IVillegas et al, 
(2010) measure systematically fainter GCLF TO for 
low luminosity galaxies. Dwarf galaxies in Virgo 
and Fornax clusters show systematically ^0.3 mag 
fainter GCLF TOs, and less populated globular clus- 
ter systems. Hence the GCLF method is not ap- 
plicable to dwarf galaxies fainter than Mb ~ —18. 
It is thus very surprising that the GCLF TO of 
the combined sample of 175 globular cluster candi- 
dates in 68 faint (My > —1 6 mag) dwarf galaxie s 
is M v<TO = -7.6 ± 0.11 mag |Georgiev et alJl2009h . 
very similar to the metal-poor GCLF TO of the MW. 
Corrections due to a different metallicity of the tar- 
get galaxy and the calibrator globular cluster sam- 
ples. 

The bimodal globular cluster color distribution 
is present in most giant galaxies. The metal- 
poor, blue, population has been shown to be old 



in essentially all well studied galaxies, while the 
red population can be either old and metal-rich 
or contain a somewhat y ounger cluster population 
(jBrodie and Straderl 120061 ) . All galaxies contain the 
blue peak in the globular cluster color distribution 
with a relatively shallow dependence of the mean 
color of that peak on the galaxy luminosity, while 
the fraction of red clusters and their mean color de- 



creases much faster for fainter galaxies (|Peng et al 
2006). The ubiquitous old metal-poor population 
can be used to minimize the corrections due to 
metallicity differ ences. This has been advocated also 
in the review by iKissler-Patid ((2000). 

For example, just selecting metal-poor cluster 
populations in the MW, M31 and other galaxies 
Di Criscienzo et al. ( 20061 ) find Mv,to values con- 
sistent to within 0.02 mag for the three samples. In 
addition as already mentioned in the previous sec- 
tion the us e of near-IR filters m inimize metallicity 



corrections ( Ashman et al. 1995^ 



(|2010h 



It should be noted here that IVillegas et al. 
use full globular cluster samples in each galaxy. Re- 
stricting the GCLFs to only the blue population 
might remove a possible bias due to different average 
metallicity of galaxies in the two clusters, but this 
correction is expected to be anyway small due to the 
near-IR z filter used in this study. 
The uncertainty in the age of the cluster sample. 

The method only works for the old globular clus- 
ters that had undergone the Hubble time of stellar 
and dynamical evolution in the potential of large el- 
liptical galaxies. The younger star cluster systems 
will display a fainter TO. This is somewhat counter- 
intuitive, as the younger clusters are brighter. How- 
ever, the fainter TO is due to GCLF peak having 
lower mass clusters and thus fainter average magni- 
tude. 

The discrepant, 0.22 mag smaller, relative dis- 
tance modulus between Fornax and Virgo clusters 
from the GCLF method could be explained by ~ 
3 Gyr younger globular cluste r systems in Fornax 



2010l ) 



cluster with respect to Virgo (jVillegas et al. 
galaxies. Such a trend in globular cluster forma- 
tion time is not unexpected based on hierarchica l 
galaxy formation simulations (|De Lucia et al . 2006), 
but again this points to a systematic dependence on 
environment, which makes the GCLF method a less 
accurate distance indicator. 



8 Summary 

The globular cluster luminosity function (GCLF) is an 
attractive distance indicator for early-type giant galax- 
ies. Through a simple functional form the characteristic 
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turnoff magnitude of the GCLF can be easily obtained 
and its intrinsic luminosity enables distance measure- 
ments up to ~ 120 Mpc, to the Coma cluster, and there- 
fore direct determination of the Hubble constant. 

It is recommended to apply the method to the metal- 
poor globular cluster sub-population and to use near-IR 
filters, as they provide least corrections for metallicity 
differences between cluster samples. The calibration 
of the V band TO magnitude for metal-poor clusters 
in the MW, M31 and 17 other nearby galaxies shows 
remarkable consistency within only 0.02 mag, and the 
weighted mean value is Mv,to = —7.66 ± 0.09, cali- 
brated whh_j^sp_ec^_to_the_JjMC distance modulus of 



18.50 (|Di Criscienzo et al.l l2006). 



The accuracy of the method is impacted however by 
a potential dependency of the standard candle on en- 
vironment as it provides a systematically shorter rela- 
tive distance, by about 0.1-0.27 mag, between the Virgo 
and Fornax clusters when compared with their relative 
distances f rom the PNLF, SBF and Cepheid distance 
indicators ( Villegas et al. 2010l) . 
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